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Abstract. We examine the problem of particle acceleration at a relativistic shocks assuming pitch- 
angle scattering and using a Hartree-Fock method to approximate the associated eigenfunctions. 
This leads to a simple transcendental equation determining the power-law index, s, given the up and 
downstream velocities. We compare our results with accurate numerical solutions obtained using 
the eigenfunction method. In addition to the power-law index this method yields the angular and 
spatial distributions upstream of the shock. 
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Particle acceleration at relativistic shocks via the first-order Fermi process has been 
the subject of a great deal of analysis since the late eighties. Initial calculations using 
semi-analytical [2., 2i,4J and Monte Carlo [1] methods based on pitch-angle scattering 
suggested a "universal" power law index for ultra-relativistic shocks of 5 = 4.22. How- 
ever, recent simulations that comp ute p article trajectories in prescribed random fields 
present a more complicated picture [iS^Isi 0,191]. Both approaches involve essentially ar- 



bitrary assumptions about the particle transport properties, and the only realistic prospect 
of making progress lies in particle-in-cell simulations. Recently, these have produced 
the first evidence that relativistic shocks can self-consistently accelerate particles [16j], 
and the results appear to favour the pitch-angle scattering approach, which motivates 
us to present here an improved approximation scheme based on [4]. A simple formula 
relating the power law index s to the shock velocities that provides a reasonable fit to 
the numerical results has already been advanced by Keshet & Waxman [|5|]. But, as these 
authors note, its derivation is incomplete. In this contribution we present a simple analyt- 
ical approximation scheme, based on the variational method familiar from Hartree-Fock 
computations of atomic eigenfunctions. The results give both the angular distribution 
and the spectal index and are in good agreement with accurate numerical evaluations. 
The particle transport equation satisfied both upstream and downstream is 

where /3 is the speed of the fluid with respect to the shock front (in which the distri- 
bution is assumed stationary) and F = (l — /3^) . We look for solutions of the form 



f{p,jl,z) = Y.ifiP ^Qiil-i) exp(A,z/r) with eigenfunction-eigenvalue pairs satisfying: 



It has previously been shown that the upstream angular distribution can be approximated 
very accurately by the eigenfunction Qi of smallest positive eigenvalue, which, in the 
ultra-relativistic case F ^ 1, is approximately Qi ~ exp( — (1 + /i)/(l — /3)) [4]. The 
value of the spectral index follows when the boundary condition far downstream is 
applied to this function, after is has been transformed into the downstream frame. This 
requires knowledge of Q\ at the downstream speed /3 = /3d, which is not ultra-relativistic. 

Using this information we chose a one parameter trial function for the upstream 
eigenfunction of the form Q\ = exp(aK/i) . a„ is determined by the requirement that Qi be 
orthogonal to the isotropic eigenfunction of zero eigenvalue, Q = constant: exp(2a„) = 
((3[,(/3 — 1) — 1)/ (a;,(/3 + 1) — 1). The corresponding eigenvalue is found by projecting 
Eq. © onto 2 1, giving 

-2al + l)+au [e^"" - l) 
^1 =^e4a„ (2a„ (/3 + 1) - 1) - (2a, (/3 - 1) - 1) ' 

which differs from the numerically calculated eigenfunction by less than 0.2% (this 
occurs near /3 = 1/3). For the downstream eigenfunction, choose a two-parameter poly- 
nomial trial function: 

Qi = l+adiPL^ad2pi'^. (4) 

As in the Hartree-Fock method, the parameters are evaluated by minimising the eigen- 
value. A minor subtlety arises here, because Eq. Q has a family of eigenfunctions of 
negative sign. However, these are easily excluded in the minimising prodedure. Enforc- 
ing the orthogonality condition and minimising the eigenvalue we obtain 

/«rfi+3j8\ 13 + 10/32 - ^25 - 20/32 + 100/34 
«.2 = -^^g— J a.i = -/3 2(-3 + 5/32) ' 



This trial function is only valid for /3 < y 3/5. The corresponding eigenvalue is 

Ai = -D ^ . ^, (6) 

lOarfi + 6adiad2 + /3 (15 + 5a^^ + 3a^2 + ^^<^d2) 

which again closely agrees with the numerically calculated value. 
The spectral index s is the solution of the integral equation [0] 

where the indices u and d refer to quantities upstream and downstream respectively. 
With the substitution t = —au{ \ + firAl^u) I Aei' the integral reduces to 



f'^ {PQ + Pit-VP2t^+Pit^) e'' dt^O (8) 

Jti 




FIGURE 1. Left: the numerically obtained eigenfunction compared to the single-parameter approxima- 
tion to Qi for P — .99. H, is the angle between the shock normal and the particle momentum measured 
in a frame of reference (effectively the downstream frame) moving at j3 = .33. Right: the numerically 
obtained eigenfunction compared to 2i ~ 1 + a^/i^u + 0^2^^ for )3 = .33. 
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(^dlAel + 2arf2. 

^2 = I {adiPiei + ^adi) n 
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^3 = ( Aei + «di Aei + adi) /al . 
Therefore, s satisfies 

P^{T{s-l,h)-T{s-\,t2))+P2{T{s-2,ti)-T{s-2,t2)) + 

i'i(r(^-3,ri)-r(^-3,?2))+/'o(r(^-4,?i)-r(^-4,f2))=o, 



(9) 



(10) 



(12) 



(13) 



where T{x^t) is the incomplete gamma function. This transcendental equation can be 
solved easily with any root finding algorithm. 

In figure [U we compare the eigenfunction Q\ obtained numerically using the Priifer 
transformation |J, 3] with our Hartree-Fock approximation for two flow speeds typical 
of upstream and downstream conditions. Agreement is very good for both /3 = .99 and 
j8 = .33. In the /3 = .33 case the approximate eigenvalue differs from the eigenvalue 
obtained numerically by about 0.2%. This inaccuracy has a small effect on the value 
of the ultra-relativistic spectral index obtained, s — 4.23, about 0.23% higher than the 
accurate result of s — 4.22, obtained using numerical evaluation of the eigenvalues and 
expanding the distribution function to higher order than 2i [|4fl. 
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FIGURE 2. The spectral index obtained for the Jiittner-Synge equation of state via the Qj method, 
Hartree-Fock approximation and the fit 5 = (3j3„ — 2/3„/3j + j3j) / (j3„ — jS^/) are compared. 

In figure [2] we compare the value of s obtained from tlie Hartree-Fock approximation 
with accurate numerical results and with the formula proposed by Keshet & Waxman 
[[st] 5 = (3/3m — 2/3;,/3j + /3j) / (/3m — jS^). In this example, the upstream and downstream 
speeds are related by the jump condition given by the Jiittner-Synge equation of state. 
The value obtained via the fitting formula is least accurate for mildly relativistic shocks 
(.5 < Fm < 2), where it deviates by Ai' ^ 0.03. The Hartree-Fock approximation deviates 
by at most A* = .015 for all shock velocities. 
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